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Abstract 

We develop a general scheme for the use of Fermi operators within 
the framework of integrable systems. This enables us to read 
off a fermionic Hamiltonian from a given solution of the Yang- 
Baxter equation and to express the corresponding L-matrix and 
the generators of symmetries in terms of Fermi operators. We 
illustrate our approach through a number of examples. Our main 
example is the algebraic Bethe ansatz solution of the Hubbard 
model in the infinite coupling limit. 
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Introduction 

The purpose of this article is to consider graded vector spaces and the graded 
Yang-Baxter algebra in such a way that the explicit construction of inte- 
grable models in terms of Fermi operators becomes easy. We will present 
simple general formulae, which will enable expression of Hamiltonian, L- 
matrix, Yang-Baxter algebra and generators of sjTiimetries in terms of Fermi 
operators, once a solution of the Yang-Baxter equation is given. 

The material developed below has its origin in two articles of Kulish and 
Sklyanin 0, where fundamental, graded integrable lattice systems were 
considered for the first time. Although most of the basic ideas are outlined in 
[0, 0, these articles remain rather sketchy as far as concrete representations 
in terms of Fermi operators are concerned. We shall try to explain in the 
following that this is a topic of its own interest. 

We would like to emphasize that the construction of fermionic represen- 
tations is a subject which has to be seen separately from the construction 
of integrable models invariant under Lie superalgebras or their deformations 
(cf e.g. 0, |, ^). Although Lie superalgebra invariant models are most nat- 
urally represented in terms of Fermi operators, there is no need to do so. On 
the other hand, there are models, such as the Hubbard model, or its infinite 
coupling limit, to be treated below, which are not invariant under graded 
algebras, but have their most interesting interpretation in terms of Fermi 
operators. The possibility of connecting a given solution of the Yang-Baxter 
equation with a fermionic representation is only restricted by a rather weak 
compatibility condition (eq. (P^D). Thus there may be different fermionic 



representations of the same model, which correspond to different gradings. 

One possibility to connect spin models with fermionic models, which is 
frequently encountered in the literature, is the Jordan- Wigner transforma- 
tion. This transformation is mostly applied on the level of Hamiltonians. 
Applying it to the Yang-Baxter algebra is a feasible, yet cumbersome task. 
The Jordan- Wigner transformation led to important progress, e.g. in the un- 
derstanding of the Hubbard model 0|. We shall argue however that it is 
an unnecessary element for the construction of fermionic representations and 
that the route taken in the present article leads to a clearer understanding of 
the issue. The Jordan- Wigner transformation does not preserve the bound- 
ary conditions. It may obstruct symmetries, and is not easily generalizable 
to an arbitrary number of internal degrees of freedom. 
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The plan of the present paper is as follows. We first present some nec- 
essary mathematical preliminaries. We introduce the parity, a concept of 
odd and even on the basis of a finite dimensional vector space, which is then 
called a graded vector space. Later this will allows us to make contact with 
Fermi operators. We extend the concept of parity to endomorphisms of the 
vector space and to tensor powers of endomorphisms. The central definition 
is equation ([T^) , which provides an embedding of a "local" basis {e^} of 
endomorphisms into a "chain of L sites", {e^} — > {e,,-^}, in such a way that 
the Cj^, depending on the value of a and (3, either commute or anticommute. 
The Cj^ are "graded analogs of spin operators". It will turn out later that 
they can be directly replaced by fermionic projection operators. This is an 
advantage of considering the grading on a basis rather than on coordinates. 
Before turning to fermions, we have to develop the quantum inverse scat- 
tering method using the basis {cj^^}. We try to keep the analogy to the 
non-graded close as possible and try to emphasize in our presentation 

the modifications due to the grading. We discuss the general case of global 
symmetries, which, due to the grading are subalgebras of Lie super algebras. 
Later, after having introduced Fermi operators, we look at the example of 
gauge transformations, where we also discuss the local case. 

Finally, we illustrate our approach with a number of examples. Most 
of these examples will be familiar to the reader. They were chosen for this 
reason. The last example, however, is new and of its own interest. Our 
approach allows us to identify the i?-matrix of the Hubbard model in the 
infinite coupling limit as a member of a family of i?-matrices of su(N) XX 
chains, which was recently proposed by Maassarani and Mathieu Since 
the example is new, we work it out in some detail. We perform an algebraic 
Bethe ansatz for the model and investigate its symmetries and the question 
of completeness of the algebraic Bethe ansatz eigenstates. 



Graded vector spaces 

Within the formalism of second quantization a physical lattice system can 
be entirely described by a set of creation and annihilation operators a^, aja 
of particles at site j of the lattice. Depending on the index a these particles 
may be bosons or fermions. Accordingly, they commute or anticommute at 
different sites, 

ajaakfn ± akfsaja = , j ^ k . (1) 
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Let us consider systems with a finite number of states per site, say m bosonic 
and n fermionic states, a = 1, . . . , m + n. We can look for matrix representa- 
tions of the operators aja acting on tensor products of the local space C"*''"". 
For spinless fermions, for instance, such a representation is provided by the 
Jordan- Wigner transformation. Unoccupied sites are bosonic, occupied sites 
are fermionic. 

If we want to describe the general situation, it turns out to be useful to 
reverse the above reasoning. Let us start from a local space of states V, 
which is then isomorphic to C™^", and let us impose an additional structure, 
the parity, from the outset. Let = Vq © Vi, and call f q G Vq even, Vi G Vi 
odd. The subspaces Vq and Vi are called the homogeneous components of V . 
The parity p is a function Vi — > Z2 defined on the homogeneous components 

of y, 

p{vi) = i , i = 0, 1 , G . (2) 

The vector space V equipped with this structure is called a graded vec- 
tor space or supcrspace. Fix a basis {ci, . . . ,em+n} of definite parity. Let 
p{a) := p{ea)- Since we want to construct an algebra of commuting and 
anticommuting operators, we have to extend the concept of parity to opera- 
tors in End(l^) and to tensor products of these operators. Let G End(l^), 
e-y — S^Ca- {ef } is a basis of End(y). If we represent as a column matrix 
with only non- vanishing entry 1 in row 7, then is an [m + n) x [m + n) 
matrix with only non-vanishing entry 1 in row a and column Define the 
parity of rows and columns, 

Pr{ei)^p{a) , Pc(e^)=p(/3) (3) 

and the parity of matrices 

p(ef)=p(a)+p(/3) . (4) 

With this definition End(y) becomes a graded vector space. 

Take X, Y from the homogeneous components of End(y), and define the 
superbracket 

[X, y]± ^XY- {-i)p(^my)yx . (5) 

Extend the superbracket linearly in both of its arguments to End(\^). Then, 
End(V^) endowed with the superbracket becomes the Lie superalgebra 
gl(m|n). Note that the above definition of a superbracket makes sense in 
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any graded algebra. We will also use it in the context of graded tensor 
powers of End(V"), which is our next issue. 

The notion of a grading may be extended to the L-fold tensor product of 
End(y), setting 

Pr{e^a\ ® ■ ■ ■ ® e^i) = + , (6) 

Pc(eS®...®e^^) = + , (7) 

p(e^; ® ...®e^^) = p{ai) + p{Pi) + . . . + p{aL) + p{Pl) ■ (8) 

It can be seen from the last formula, that homogeneous elements A = 
^%..1^t^ai ® ■ ■ ■ ® of (End(V))®^ with parity p{A) are characterized 
by the equation 

(_l)E^=i(pK)+f(/5.))^«i---«^ ^ . (9) 

This implies that AB is homogeneous with parity 

p{AB)^p{A)+p{B) , (10) 

if A and B are homogeneous. 

The above definitions allow us to introduce another tensor product struc- 
ture ®s, which is called the graded or supertensor product, on exterior 
powers of End(V"). Choose v and w from the homogeneous components 
of (End(y))®'^' and (End(y))^', respectively. Then, by definition, 

^; ®^ ^Z; = ^; ^i; . (11) 

As a simple consequence of this definition, the supertensor product is asso- 
ciative, {u <^s v) <^s w — u <S>s {v <S>s w). This follows first for u, v, w taken 
from the homogeneous components of certain powers of End(y), and then 
by linearity for arbitrary u, v and w. 

Here comes the central definition of the paper. The supertensor product 
induces an embedding of into (End(V))'^^, 

p./3 _ r^-O-l) (5?) p/3(5?) T®s{L-j) 

^ ^_^^(pia)+pm EL,+if(^'=)/^5-^i) ® ® e:^^+^ ® . . . ® e^^ . (13) 
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Im+n in this equation denotes the (m + n) x (m + n) unit matrix. In the 
second equation summation over double tensor indices is understood. The 
index j on the left hand side of (|T2]) will be called site index. The matrices 
Cj'^ realize relations of the form (|l]). For j k we find 

eife4 = (-l)(^(")+^(/^))(^W+p(^))e4e,f . (14) 

It follows from eq. ([131) that e^^ is homogeneous, and that 

pie/j = pia) + pi/3) . (15) 

Hence, in agreement with intuition, eq. ( p!^ says that odd matrices mutually 
anticommute, whereas even matrices commute with each other as well as 
with the odd matrices. For products of matrices e^^ which are acting on the 
same site (|13]) implies 

e/e/ = 5^e/ . (16) 
Using the superbracket, (|T^) and (|1^) may be combined to 




The right hand side of the latter equation with j = k gives the structure 
constants of the Lie superalgebra gl(m|n) with respect to the basis {ej(^}. 



The permutation operator 

The permutation operator plays an important role in the construction of 
local integrable lattice models. In the graded case it requires the following 
modifications of signs, 

P,k = {-ir^'^e/^e,1 . (18) 

This operator induces the action of the symmetric group 5*^ on the site 
indices of the matrices ej^. The following properties are easily verified, they 
follow from (|T|) and (|16D" 

[m — n)id , 
-- id , j^k , (19) 



(a) 


Pkj 


(b) 


P 


(c) 




id) 


Pjk^ka 
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Pjkcd 
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eiiP.u , Ji^l^k 
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Figure 1: The Yang-Baxter equation is most easily memorized in graphical 
form 

Because of {d) and (e) Pjk generates a faithful representation of 5"^, 



From the right-hand side of this equation we can read off the matrix elements 
of Pi2 with respect to the canonical basis of End(y ® V). 

The graded Yang-Baxter algebra 

In the present context it is most suitable to interpret the Yang-Baxter equa- 
tion as a set of functional equations for the matrix elements of an (m-|-n)^ x 
(m + n)^-matrix R{u, v). We may represent it in graphical form as shown in 
Figure 1, where each vertex corresponds to a factor in the equation 



(m, v)R%7^, {u, w)Rp^n {v, w) = Rp^, {v, w)R'^^n {u, w)R'^fp„ {u, v) . 

(22) 



Note that there is a direction assigned to every line in figure 1, which is 
indicated by the tips of the arrows. Therefore every vertex has an orientation, 
and vertices and /^-matrices can be identified according to figure 2, where 
indices have been supplied to a vertex. Summation is over all inner lines in 
figure 1. 




(20) 



Let L = 2. Then 
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Figure 2: Identification of tlie i?-matrix with a vertex 



Starting from the Yang-Baxter equation we will construct a fundamental 
graded representation of the Yang-Baxter algebra. For comparison let us 
briefly recall the non-graded case = 0). Define R{u,v) as 

R^Jiu,v) = R'^nu,v) . (23) 

Introduce the L-matrix at site j, 

ir.M = RjlMe.t, . (24) 

Then multiplication of the Yang-Baxter equation ([22| ) by Cj^" implies that 

R{u,v) {Lj{u,w) ® Lj{v,w)) = {Lj{v,w) ® Lj{u,w)) R{u,v) , (25) 

where the tensor product is now a tensor product between matrices, according 
to the convention {A (g) B)'^] = A'^B]. We may replace Lj in (|25D by some 
matrix T and may interpret it as defining an abstract algebra for the matrix 
elements of T. This algebra is called a Yang-Baxter algebra with i?-matrix 
R. Lj is called its fundamental representation. Since L-matrices at different 
sites commute, any product of L-matrices with different site indices is another 
representation of the same Yang-Baxter algebra. 

The construction of a graded Yang-Baxter algebra and its fundamental 
representation requires only minimal modifications of the above scheme. Let 
us assume we are given a solution of (0), which is compatible with the 
grading in the sense that 

i?°f(M,t;) = (-l)P(-)+f(^)+fW+P(^)i?°f(M,t;) . (26) 
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Define a graded L-matrix at site j as 

/:,^(«,t;) = (-ir(")^Wi?^J(«,T;)e4 . (27) 

Eq. ( p6D implies that the matrix elements of Cj{u,v) are of definite parity, 

p{C,';{u,v))=p{a)+p{(3) , (28) 

and that they commute as 

C,';{u,v)Ck]{w,z) = (-l)(^(")+^('3))(PW+^(^))£j(t^,^)£,;(w,t;) . (29) 

It further follows from the Yang-Baxter equation (^21) and from eq. ( p6| ) that 

R{u,v) {Cj{u,w) Cj{v,w)) = {Cj{v,w) <^s '^j{u,w)) R{u,v) . (30) 

In analogy to the non-graded case above, the supertensor product in this 
equation is to be understood as a supertensor product of matrices with non- 
commuting entries, {A ®, 5)^J = 1-1)(p('^)+pW)pM A'^B] . In a sense this 
definition is a contravariant counterpart of equation (|11]). Given matrices A, 
B, C, D with operator valued entries, which mutually commute according 
to the same rule as Cj and Ck in eq. (P^, we obtain for the product of two 
supertensor products 

{A ®s B) {C ®s D) = AC ®s BD . (31) 

Eq. (|30|) may be interpreted as defining a graded Yang-Baxter algebra with 
i?-matrix R. Cj is then its fundamental representation. 

Starting from (|30| ) we can construct integrable lattice models as in the 
non-graded case. Let us briefly recall the construction with emphasis on the 
modifications that appear due to the grading. Define a monodromy matrix 
T(m, v) as an L-fold ordered product of fundamental L-matrices, 

T{u,v) = Cl{u,v) . . .Ci{u,v) . (32) 

Due to eq. ( p!0D the matrix elements of T(m, v) are homogeneous with parity 
p{T^{u,v)) = p{a) + p{P). Repeated apphcation of (|30|) and (^) shows 
that this monodromy matrix is a representation of the graded Yang-Baxter 
algebra, 

R{u,v){T{u,w)(g)sT{v,w)) = {T{v,w)^sT{u,w))R{u,v) . (33) 
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In the non-graded case (n = 0) the supertensor product in p3| ) agrees with 
the usual tensor product. Multiplying (|3^) by R~^{u,v) and taking the 
trace of the whole equation then implies that [tT(T{u,w)),tT(T{v,w))] = 0, 
and the transfer matrix tT{T{u,w)) provides a generating function of mu- 
tually commuting operators, which may take the role of conserved quanti- 
ties of an integrable lattice model. For non-trivial grading the trace has 
to be replaced by the supertrace, which is generally defined as str(74) = 
(_l)E,liPK)^m...j^iv_ Then (|3|) implies that 



[str(T (m, w)), str(T {v, w))] = , (34) 

in complete analogy with the non-graded case. 

Let us assume that R{u, v) is a regular solution of the Yang-Baxter equa- 
tion. This means that there are values Mq, "^0 of the spectral parameters such 
that R?;^iuo,Vo) = 626^. Then (P) implies 

^■;K,t;o) = (-l)^(")^(^)e,^ , (35) 

and we can easily see that the supertrace of the monodromy matrix at {uq, fo) 
generates a shift by one site, 

str(r(«o,^o)) = {-ir^^K"{uo,vo) 

= Pi2P23 ■ ■ ■ Pl-il =■■ U . (36) 

This implies that t{u) := ln(str(T('u, tig))) generates a sequence of local 
operators which, as a consequence of (|3^) , mutually commute, 

t{u) =ifl+{u-uo)H + 0{{u-uof) . (37) 

n in this expansion is the momentum operator. On a lattice, where the 
minimal possible shift is by one site, and thus U rather than 11 is the fun- 
damental geometrical operator, some care is required in the definition of 
n. As was shown in a proper definition may be obtained by setting 
n := — iln(f/) mod27r and expressing the function f{x) = a;mod27r by its 
Fourier sum. Then 11 becomes a polynomial in II. 




^ = ^T. 7; + -n^] , (38) 
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where = l-njL. The first order term H in the expansion ( p7| ) may be 
interpreted as Hamiltonian. Using (|36|) it is obtained as 



ii+i 



where i^LL+i = i^Li and 



(39) 



(40) 



We would hke to draw the reader's attention to the following points, (i) 
The i?-matrix R in equation ( pOD does not undergo a modification due to 
the grading, (ii) The only necessary compatibility condition which has to be 
satisfied in order to introduce graded L-matrices is equation (^). As we will 
see in the examples below, this is a weak condition. A given i?-matrix may 
be compatible with different gradings, leading to different graded L-matrices. 



A first example 

As mentioned above the matrices e^^ for fixed j form "local" representations 
of gl(m|n). After summing over all sites, we obtain a "global" representation, 

K = Ee.f , (41) 
[^f,E^]± = - (-1)(pH+p(/3))W7)+p(-5)) ^5^^ _ (42) 

The are symmetric in the site indices by construction. Thus [-Pj^+i, E^] = 
0, and we obtain a gl(m|n) invariant Hamiltonian, if we are able to find a 
solution of the Yang-Baxter equation (^21) , which leads to Hjj^i = Pjj+i 
in equation (EDI). Then, comparing (O) and (ffOl), duR'^s(u,vo) = 

' ' — ' ' — ' ' u=uo 

(^—iyi°')pWS'^6^. Taking into account regularity, we find the following mini- 
mal Ansatz for R{u,v), 

Rf&{u, v) = 5^5^^ + {u- v){-lY^'^^^^P^5^5l , (43) 
which is indeed a well known rational solution of the Yang-Baxter equation 
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Global symmetries from local symmetries 

We are going to consider now the general case of symmetries of the mono- 
dromy matrix, which stem from Lie superalgebra invariance of the i?-Matrix 
in a sense to be specified. Choose homogeneous from gl(m|n). Let 

xj ■= x^e/^ and X := E}=iXj. Define R"p]{u,v) := v). 
Assume that R{u, v) satisfies the invariance equation 

R'^7s{u, v) x^p + FQ,{u, v) xf 

From here we can move step by step to the invariance of the transfer matrix 
str(T('u, f )). Contraction of ( PD with e^^ yields 

and it can be shown by induction over L that the monodromy matrix satisfies 

T^^{u,v)x} + T^{u,v)X = 

(_l)P(-)(p{7)+P(/3)) a;°r/(M, v) + (_i)p(-)(p(°)+p(/5)) XT^{u, v) . (46) 

We finally take the supertrace of this equation and arrive at 

[str(T(M,t;)),X] = . (47) 

In many cases the symmetry of the i?-matrix is evident by construc- 
tion, e.g. when the i?-matrix is an intertwiner of representations of quantum 
groups. Yet there are examples, as Shastry's i?-matrix of the Hubbard model 
[^, 1^, where the symmetries are less obvious [jlOl , Moreover, as can be 
seen by the above derivation, the symmetries of the transfer matrix are deter- 
mined by the symmetries of R rather than R. The symmetries of R depend 
on the choice of the grading. 

It may be argued that, in presence of a grading, the matrix R is more fun- 
damental than i?, since R determines the L- matrix, eq. (|2^) , the symmetries 
of the model and, if it exists, the semi-classical limit pi]. Substituting R into 
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the Yang-Baxter equation (|22|), we obtain the so-called graded Yang-Baxter 
equation, which equivalently might have been taken as the starting point of 
our section on the graded Yang-Baxter algebra. Since it is the non-graded 
matrix R, however, which fixes the structure of the Yang-Baxter algebra, eq. 
(|30D, we stood away from this point of view. 



Representations in terms of fermions 

In this section we shall explain how the various graded objects, which have 
been introduced so far, can be expressed in terms of Fermi operators. To 
begin with, consider spinless fermions on a ring of L lattice sites, 

{cj.Ck} = {c],c\} = Q , {cj,c\} = 5jk , j, A; = 1,...,L . (48) 

Locally there are two states, every site is either occupied by a fermion or 
it is empty. Slightly deviating from the usual language we may say that Cj 
and c] annihilate or create the occupied state. Let us define a pair a^-^, aj\ 
of (trivial) annihilation and creation operators of the unoccupied state by 
setting aj^ = ajl = 1, and let us write aj^ = Cj, ajl = c^j- Let rij = c^jCj 
denote the density operator. The operators 

Y/^ = a,i{l-n,)a,p (49) 

are then obviously local projection operators, i.e. they satisfy 

y /3y / = s^Y/ . (50) 

They carry parity, induced by the anticommutation rule for the Fermi oper- 
ators. Let j 7^ k. Then 1^-^ and Y^^ anticommute, if both are build up of an 
odd number of Fermi operators, and commute in all other cases. This fact 
can be expressed as follows. Let p{l) = 0, p{2) = 1 and p{Yj^) = p{a) +p(/3). 
Then Yj^ is odd (contains an odd number of Fermi operators), if p{Yj^) = 1, 
and even, if p{Yj^) = 0. The commutation rules for the projectors Yj^ are 
thus 

Yj%'^ = (-lYp(»)+PimpM+PiS))Y,^^Y/^ , (51) 

As a memorizing scheme for the projection operators let us combine them 
into the matrix (Yj)'^ = Yj^, 

( , 



\ 4 



(52) 



14 



0]) and (|51| ) are representations of equations (|T3|) and ( ]T6| ) in the particular 



case m = n = 1. Since all our considerations in the previous sections en- 
tirely relied on eqs. ( [T^ ) and (|TBp, we may simply replace e^^ by Y^-^ in eqs. 

and (|40|) to obtain fermionic representations of the permutation 



operator, the L-matrix and the Hamiltonian. The permutation operator (for 
j 7^ k) becomes 

p V ^ J- V 2 V 2 V ^ V 2 V 2 

= l-(ct-4)(c,-c,) . (53) 

Fermionic representations compatible with arbitrary grading can be con- 
structed by considering several species of fermions and graded products of 
projection operators. We shall explain this for the case of two species first. 
This is the most interesting case for applications, since we may interpret the 
two species as up- and down-spin electrons. We have to attach a spin index 
to the Fermi operators, Cj — > Cjo-, a =],i, {cja,cl^} = ^jk^ar- Accordingly, 
there are two species of projection operators, Yj^ — ^ ' 

yT^yi"^ = ('_i)(p{")+p(/3))(p(7)+P(5))yiYT^ _ 

Let us define projection operators for electrons by tlie tensor products 

Y/J, = (-l)(^H+^(/^)M^)l.^V;; = {yI ®, Yjyj^ . (55) 

Then 

y^-^^ inherits the parity from Y^^^ and Yj ^- The number of Fermi operators 
contained in Y/^ is the sum of the number of Fermi operators in Y}^ and 
F/^. Hence piYj^J^) = p(r/^) + p(r/^) = p{a) + ...+ p{6), and the analog 
of holds for Y/J^, too. A gain we present all projection operators in form 
ofamatrix(y,)^J = y,f^, 



15 



%i)(l-^jT 
"-it) 



V 



I c 



.(57) 



Here we used the standard ordering of matrix elements of tensor products, 
corresponding to a renumbering (11) 1, (12) 2, (21) — > 3, (22) — > 4. 
Within this convention Y^|^!^ is replaced by Yj^, a, (3 = 1, ... ,4, which then 
satisfies (pOD and (0) with grading p(l) = p(4) = 0, p{2) = p{3) = 1. 

The permutation operator of electrons turns out to be a product of per- 
mutation operators of up-and down-spin electrons, 



jk 



i-)p(/3)+p('5)y./3'5y "7 

.^^p(/3)+p(^)+(p(a)+p{/3)){p(7)+p(5))yi/3yT'5yi°yTT 



-1) 



pj. pT 



pT pi 



(58) 



So far we have considered the case of spinless fermions with two-dimen- 
sional local space of states and grading m = n = 1, and the case of electrons 
with four-dimensional space of states and grading m = n = 2. There are four 
different possibilities to realize ([l^ ) and (|l^) in case of a three-dimensional 
local space of states, m + n = 3. They can be obtained by deleting the a's 

1,2,3,4. dO) and (||) 



a 



row and column of the matrix Yj in eq. { ^7\ 
remain valid, since the operators Yj^ are projectors. 

As an example let us consider the case where the fourth row and column 
of Yj, eq. (^7|), are deleted. The local Hilbert space is then spanned by the 
three states |0), c]||0), c]jO), i.e. double occupancy is now excluded. The 
operator 

E = 1 - r> = 1 - n,^n,^ (59) 

projects the local Hilbert space of electrons onto the space with no double 
occupancy. The global projection operator for a chain of L sites is given by 
the product 



(60) 
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The permutation operator Pjk is again given by eq. ( |I8| ) with Yj'^ replacing 
Cj^ . Summation is now over three values, a,(3 = 1,2,3, and the grading 
is p(l) = 0, p{2) = p{3) = 1. An elegant way of taking into account the 
simplifications arising from the restriction to the Hilbert space with with no 
double occupancy is to consider PjkA instead of Pjk- Since rij-^njiA = 0, we 
obtain 

PjkA = A{clck^ + cLc,.) A - 2{S;St - \n,nk)A + (1 - n, - n^) A . (61) 

Here we have introduced the electron density rij = rij^ + n.ji and the spin 
densities 

S'; = \<pcy^P ■ (62) 
The a", a = x,y, z, are the Pauli matrices, and we identify 1 with | an 2 
with I in the summation over a and j3. The spin densities can alternatively 
be written as 

S] = H^%Y/:;l ■ (63) 

Using 

{a'^^a%] = 2626}-6-,6] , (64) 

we obtain 



a,P=2 j 

(65) 

which gives the second and third term on the right hand side of (|6TD. Note 
that the permutation operator Pj^ in eq. (|6TD is no longer a product like in 
eq. dH). 

It should be clear by now, how to generalize the above considerations to 
an arbitrary number of species of fermions. In the case of iV species we may 
define 

Then 

y/3i.../3^Y-5i...5^ = , (67) 

Jai...ojv -'71---7JV 71 7iV Ja\...aiq ' V / 

■i^/3i.../3iv V 5i...5iv _ / 1 \Z]ffc=i(PK)+P(/^j))(P('>fc)+P(^fc))v;<5i---5iv iao\ 
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which can be shown by induction over the number of species. Here the 
grading is m = n = 2^~^. The most general case is obtained by deleting 
rows and columns from Yj, eq. (^Bj), in analogy to the example considered 
above. 

Let us note that the operators Yj^ for two species of fermions appear 
under the name Hubbard projection operators in the literature. 



Gauge transformations 

The canonical anticommutation relations (^Sp are invariant under local gauge 
transformations 



Zj , e [0, 27r] . (69) 



Thus all our formulae remain correct, if we replace Cj by Cj. The local 
gauge transformation induces a transformation of the matrix Yj of projection 
operators, 

Yj — > Y^- = e''''^'°"^/^Fje~''^^''^^''^ = e~^'^^"^y^-e''^^'"^' . (70) 

Here is a Pauli matrix. In the case of N species of fermions there are 
gauge parameters y^J, a = 1, . . . ,N . Since e^"^^" ^"^ is diagonal, the form ([70|) 
of the transformation rule for Yj carries over to the case of species. Let 

G, = e<"'/2 ^ . . . ® e'^^'/2 ^ (71) 

Then 

r, ^ Y, = G,Y,G] . (72) 

For the sake of simplicity we assume in the following that we are dealing with 
two species of fermions. Note, however, that the following considerations also 
apply in the most general case, where we would have to deal with a sub-matrix 
of Yj, eq. (|66D. The L-matrix (|27|) behaves under gauge transformations as 

A;=G]:,(-i)'"°'""-''fii°:;'(«.^)i'4'G.? . (73) 

Here we defined 

Rj{u, v) = (G, ® G,)R{u, v) [G] ® G]) . (74) 
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Let us consider two implications of gauge invariance of the i?-matrix. 

(i) Global gauge transformations: Assume that Gj = G, j = 1, . . . L, and 
that R{u,v) = R{u,v), say, for arbitrary ip\ ip^ (|= 1, |= 2). Taking the 
derivative of (|74D with respect to at ip'^ = ip^ = yields 

[R{u,v),l2(^(^^®h + h^h(^(r"] = . (75) 

Since /2 cr^ /4 and 1^® I2® cr^ are diagonal, we may replace R{u, v) in 
this equation by R{u,v). Then (^ip is satisfied with x = I2® cr^. Thus the 
transfer matrix commutes with 

X = jz{Y,\-Y,l + Y,l-Y,f)=j2{l-2n,,) , (76) 
i=i i=i 

which means that the number of up-spin electrons is conserved. A similar 
statement is easily verified for the number of down-spin electrons or, in the 
weaker case, when (|7^ is only satisfied for Lp^ = for the total number of 
electrons. 

(ii) Local gauge transformations: Assume that Gj = {GY , where G is a 
fixed matrix as in the example above, and that R{u, v) is invariant under G. 
Then 

C,{u,v) = {G^y Cj{u,v){Gy , (77) 

and the transformed monodromy matrix becomes a simple expression in 
terms of the original L-matrices, 

f{u, v) = {G^Y Cl{u, v)GCl-i{u, v)G... Ci{u, v)G . (78) 

Transformations of this type can be used to introduce a phase factor into a 
typical nearest-neighbor hopping term as it appears in the Hubbard Hamilto- 
nian. If the phase factor is just —1, then the hopping term changes sign. The 
factor (G^)^ generally modifies the transfer matrix str(T('u, v)). It leads to a 
twist of the periodic boundary conditions. Note however, that it may happen 
for certain values of L, e.g. if L is divisible by 2 or by 4, that (G^)^ = 1. 

Gauge transformations of the type considered modify the shift operator. 
Using the right hand side of ( [701 ) in eq. ( PB| ) we can easily see that 
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where summation over a =t, i is implied. The first term in the exponent of 
(|79D is related to a twist of boundary conditions. 

At this point will will not go into further detail, because we think that 
more detailed considerations are only sensible in the context of concrete 
models. Let us only remark that global symmetries of a given model may be 
affected by the exponent in eq. ([79|), even if there is no twist, i.e. even if the 
first term is equal to zero modulo 27r. 

In the context of integrable fermionic models the term gauge transfor- 
mation has, unfortunately, two different meanings, which should not be con- 
fused. Besides the meaning discussed above it is also used for transformations 
on the i?-matrix of the form R ^ {G ® G)R{G~^ (8) where G is not 

necessarily a diagonal matrix. Transformations of this form leave the Yang- 



Baxter equation (p2D invariant. 



More examples 

To illustrate the formal discussions of the preceding sections let us further 
elaborate on the models with rational /^-matrices (|43|) . These models have 
been studied exhaustively in the literature. They appeared first as lattice 



gas models of Lai [T^ and Sutherland [0 and were solved by coordinate 
Bethe ansatz. Later Kulish studied them within the framework of the graded 
quantum inverse scattering method and obtained spectrum and eigenstates 
by means of the nested algebraic Bethe ansatz Still Kulish did not write 
down any Hamiltonian density in terms of Fermi operators. This was first 



accomplished by Schlottmann for the gl(l|2) invariant case |]T4|. The Hamil- 
tonian in fermionic representation is the Hamiltonian of the super symmetric 
t-J-model. Schlottmann solved it again by means of the coordinate Bethe 
ansatz. However, he was not aware of the underlying algebraic structure. 
The underlying algebraic structure was successively unravelled by different 
authors |T3|, |TB|, leading to a solution of the super symmetric t- J-model by 



nested algebraic Bethe ansatz |1T7| , |T8[. The gl(2|2) invariant model with the 



Hamiltonian in fermionic representation was studied by EBler et al. []T^ pp] ]. 

Let us write down Hamiltonian, i?-matrix, L-matrix and generators of 
symmetries for the gl(l|l) invariant case using the fermionic representation 
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2|). The i?-matrix follows from (|43| ) with p{l) = 0, p{2) 



I 



R{u, v) 



1 + u — V 



1 u — V 
u — V 1 



V 



1 — U + V J 



The L-matrix is obtained from (12 

/ 



Cj{u,v) 



u — V + e 



\ 



^J2 



u — V — e 



3 2 



U — V + 1 ~ Tlj 



U — V — Tin 



As Hamiltonian density we may choose 



-P. 



jj+l — {Cj - Cjj^i){Cj - Cj+i) - 1 

- 1 



(c]cj+i + c]+iCj) + rij + rij+i 



(80) 



(81) 



(82) 



(83) 



(cf (^31)), which is the Hamiltonian density of a system of free fermions. The 
odd generators of symmetries are 



El 



El 



L 



L 

E^i • 



(84) 



(85) 



From the even generators we can construct only one non-trivial combination, 



UEl~El + L) 



i=i 



(86) 



Of course, being a model of free fermions, (^) is trivial in the sense that the 
corresponding Hamiltonian can be diagonalized by Fourier transform. We 
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wrote down the above formulae in order to illustrate the general formalism. 
The reader can easily repeat all the steps for the case of the super symmetric 
t-J-model, where the Hamiltonian density is given by -f/jj+i = —Pjj+i with 
Pjj+i according to eq. (|6lD. 

Let us indicate the Hamiltonian density of the gl(2|2) invariant model 



1^, |2^ using the fermionic representation (|57D. p3|) and ( ^OD imply that 



Here the minus sign on the right hand side has been chosen in order to meet 
the usual conventions. If we open the brackets in (|87|), we obtain the rather 
voluminous expression of EBler, Korepin and Schoutens (cf e.g. eq. 

(8) in [|1^), which we do not repeat here due to space limitations. Since 
[Pjj^i,E^] = 0, a, P = 1,2,3,4, we may add, for instance, UEf, U real, 
to the Hamiltonian obtainable from ( |S7D without spoiling integrability (see 



1^, ^). This amounts to adding a Hubbard interaction Urij^rij^ to the 
Hamiltonian density. Note that the factorized form (^TD of -f^jj+i appears to 
be new. 

Here is a more recent example. Maassarani and Mathieu constructed an 
"su(A^) analogue" of the XX-model and found the corresponding i?-matrix 
[§]. We will show that a certain fermionic representation of the "su(3) ver- 
sion" of this i?-matrix generates the Hamiltonian of the Hubbard model at 
infinite coupling. Let us write the i?-matrix of Maassarani and Mathieu in 
the form 

N 



a=2 

N ( 

+ E(e}®e^ + e"®e})+cos(.-.) e}®ei+ E 



/9 



a=2 \ a, 13=2 



Note that R depends on an additional free parameter 6. In order to meet our 
conventions we modified the original i?-matrix of Maassarani and Mathieu by 
a gauge transformation G,R^ {G0G)R{G-^0G-^), Ge^G-^ = e[^'!;_i)ZtN- 
Clearly, R is regular. Let us now restrict ourselves to the case = 3. Then 
R is compatible with the grading p(l) = 0, p(2) = p{3) = 1. Let us consider 
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the fermionic representation which is obtained from Yj, eq. (|57|), by deleting 
the fourth row and column. Using (^0]) we obtain the Hamiltonian density 

Hjj+i= ^ (e''^c]^Cj+i,^ + e"''^c]+i^^Cj^)(l-nj-_^)(l-nj+i,_^) . (89) 

We see that 6 is connected to local gauge transformations. Setting 6 = n, 
Hjj^i turns into the Hamiltonian density of the Hubbard model at infinite 
coupling, which is the same as the restricted hopping part of the super sym- 
metric t — J-Hamiltonian. Hence the Hamiltonian can be written in the more 
familiar form 

L 

i7 = -^A(4c,+i,. + c]+i,,c,.)A , (90) 
i=i 

where A, eq. (^), is the projection operator which excludes double occu- 
pancy. Of course, (^) only makes sense, if the number of particles is less 
than L. 



Algebraic Bethe ansatz for the Hubbard model in the 
infinite coupling limit 



We can now apply the algebraic Bethe ansatz ||2T| to construct eigenvectors 
of the Hamiltonian (pOl). We shall regard equation (|33|) as a set of algebraic 
relations between the elements of the monodromy matrix T{u,w). Since 
T{u, w) is a function of m — w in this case, we shall simply write T{u, w) = 
T{u — w) hereafter. Following Maassarani and Mathieu , we shall use the 
notation 

CAu) 



T{u) 



Siu) 



C2{U) 



\ 



til 



u 



tl{u) 



(91) 



tl{u) I 



Let |0) denote the vacuum state defined by Cj|0) 
of the L-matrix, which follows from (l2l 



0. Using the explicit form 



/ 



Cj{u) 



y.2 



'^3 2 



- sin uYj^ —cos uYj^ 



- cos uYj^ 



^3 3 



- cos « y^ 3 



- sin u Yj ^ —cos u Yj g 



(92) 
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we can easily see that 



THIO) 



s{u)\0) * 
t{u)\0) 








(93) 



{u)\0) ) 



where s{u) = cos-^(u), t{u) = (—1)^ sin'^(M). This equation indicates that the 
vacuum |0) is an eigenstate of the transfer matrix. Following the procedure 
of the nested algebraic Bethe ansatz p|, O, we assume the state 



\X,,---,Xp)=F^^-^-Ca,{Xl)...Ca^{Xpm 



(94) 



to be an eigenstate of the transfer matrix. Here summation over a-,- = 1, 2 is 
understood. 

The commutation rules between the elements of the monodromy matrix 
can be extracted from (ETi, 



S(u)Ca{v) = COt('U — v)Ca{v)S{u) — 



tl{u)Cb{v) = COt{u - v)Ca{v)ttiu) 



Ca{u)S{v) 

sm{u — v) 

Ca{u)tl{v) 



(95) 
(96) 

(97) 



sin('u — v) 

Comparing these commutation rules with the ungraded case in we notice 
that there appear extra minus signs when commuting ^^{u) and Ch{v). We 
use the above relations to calculate the action of the transfer matrix str(T(M)) 
on the state |Ai, . . . , Xp). The resulting terms are classified as wanted terms 
and unwanted terms, respectively. The wanted terms are 
p 

Jl cot(M - Xj){s{u)\Xi, . . . , Ap) 

-t{u){r^P^F)"^CaA^,)...CaMpm} ■ (98) 

Here r(p) is the shift operator on a 2 dimensional auxiliary space, which may 
be realized as the space of states of a p-site spin one-half chain. Its matrix 
elements are given as 



ai ...ap 



dp ai 



(99) 
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The wanted terms are proportional to |Ai, . . . , Ap), if F is an eigenvector of 
the shift operator r^^^ with eigenvalue A*^^-*. On the other hand, in oder for 
the unwanted terms to cancel each other, we have the condition 



r(p)F= (-l)^cot^(Aj)F 



(100) 



The spectrum of r(p) is fixed by the condition [t'^p^Y = 1. This condition 
implies that the spectrum consists of powers of the p-th root of unity (cf 
appendix B of The eigenvalue A(m) of the transfer matrix str(T(M)) is 

thus given as 



A{u) 



SIM — t{u)e p 



Y[ COt(M 



A,) 



where Xj, j = 1, . . . ,p, is a solution of the Bethe ansatz equations 



-l)^cot^(Aj 



i27ri 

e p 



I = 0, 



(101) 



(102) 



The only remaining problem is the explicit construction of eigenvectors of 
the shift operator t^p\ This is, of course, a highly degenerate problem. We 
may use any integrable spin one-half chain with regular i?-matrix to solve 
this problem. A particularly simple choice is the XX spin chain. Its i?-matrix 
R is given by eq. (|88| ) with N = 2 and 6 = 0. The L-matrix in (non-graded) 
spin representation is 



Lj{u) 



cos{u)ejl + sin(u)e 



2 



--J 2 



sm{u)ej\ + cos{u)ej 



\ 



2 / 



(103) 



Let us use the notation 



Tiu) 




A{u) 
\C{u) 

for the monodromy matrix T{u) = Lp{u) . . .Li{u). With the choice |0) 
^ of auxiliary vacuum, we obtain 



(104) 



T{um 



cosP(u)|0) * 

smP{u)\0) 



(105) 
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which shows that the auxihary vacuum is an eigenstate of the transfer matrix 
tr(T(M)). To construct the eigenstates of tr(T(M)) we need the commutation 
relations 

[B{u),B{v)] = , (106) 
A{u)B{v) = -cot{u-v)B{v)A{u) + ^^j^^^ , (107) 

Diu)Biv) = cot(u-v)B(v)D(u)-^^^$^^ , (108) 

smlM — v) 



which are part of the Yang-Baxter algebra (P3|). Proceeding as above, we 
consider the action of the transfer matrix tr(T(u)) = A{u) + D{u) on the 
state 

|//i,...,/i,) = 5(/ii)...5(/i,)|0) . (109) 
It turns out that . . . , fig) is an eigenstate of tr(T(M)) with eigenvalue 

s 

A(p)(u) = {(-I)^cosP(m) + smP{u)} [] cot(M - /i^) , (110) 

if the parameters /ij, j = 1, . . . , s, satisfy the equation 

cot^(/i,) = (-1)^+^ . (Ill) 
Eq. ( |110| ) imphes in particular that 

A(^)(0) = ncot(/x,) . (112) 

In order to introduce a convenient basis let us use the operators ej\ ej\ 
applied to the auxiliary vacuum |0) places an |-spin at site j, if a = 1, 
and a |-spin, if a = 2. The states ei\^ . . . e^^ |0) form an orthonormal basis 
of our auxiliary Hilbert space, since = S^ejl and ej\\0) = |0). The 

components of the shift operator U = tr(T(0)) with respects to this basis are 

(0|e//...ei^f/eii^...e,^jO) = 

(0|e,^ . . . e,Wa, . . . e'e.ljO) = 6l%\ . . . , (113) 
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such that we can identify r*^^^ with JJ . Let F = . . . ,fis)- It follows from 
eq. ([TT|) that 

s 

^{p)p ^ Yl cot{fij)F . (114) 
i=i 



Now eq. ( |1 1 Ij ) implies that 



p 



ncot(/.,) = (-1)^(^+1) = 1 , (115) 

and we have verified that the eigenvalues of r^^^ are indeed powers of the pth 
root of unity. Note the the components of F with respect to our basis, which 
enter the definition (0) of |Ai, . . . , Ap) can be written as 

(0|epi^..ei?i5(/ii)...5(/i,)|0) . (116) 



-iai...ap 



We can use the symmetries of the monodromy matrices T{u) and T{u) 
to deduce restrictions on the numbers p and s. Eq. ( PD implies that the 
i?-matrix of the XX chain is invariant under a^, [R, ® I2 + I2 ® cr^] = 0. 
Let 

S^ = E^; . (117) 
i=i 

This is twice the operator of the z-component of the total spin. We infer 
from (^6[) that [r{u),a^] = [S^,T(-u)]. B{u), in particular, commutes with 

as 

S"5(n) = 5(m)(S^ - 2) . (118) 

Thus 

|S^|/ii,...,/i,) = i(p-2s)|/ii,...,/i,) , (119) 

l/ii, . . . , /is) is an eigenvector of the z-component of the operator of total 
spin with eigenvalue |(p — 2s). This eigenvalue cannot be smaller than — |p. 
Therefore s < p. 

Using once more eq. (^^ we find that the i?-matrix of the Hubbard model 
in the infinite coupling limit is invariant under the even generators of gl(l|2), 
e\, e|, 63, 63, 63. They span the Lie algebra gl(l)©gl(2). The matrix + 
generates the symmetry operator 

^ = E (^.2 + ^.3) = EKt + nn - 2n,Tn,i) . (120) 
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This is the particle number operator on our restricted Hilbert space, where 
double occupancy of sites is forbidden. Eq. (^) implies that [A^, T{u)] = 
[T{u), el + 63] and, in particular, that 

NCa{u) = Ca{u){N + l) . (121) 

We conclude that 

iV|Ai,...,Ap) =p|Ai,...,Ap) . (122) 

The state |Ai, . . . , Ap) is constructed by acting on the vacuum with operators, 
which leave the Hilbert space with no double occupancy invariant. Therefore 
it can contain at most L particles, and p < L. 

Let us perform a naive counting of states. We will start with the auxiliary 
spin states . . . , fig)- Eq. ( |lllj ) has p solutions for every j = 1, . . . s. The 
equivalence of XX chain and a chain of spinless fermions suggests that all 
fij have to be different. States . . . , /i^), which merely differ by the order 
of the Hj, are identical due to ( p,06| ). It follows that we can order the fij as 
Hi < . . . < fig- The number of solutions of ( |111| ) which satisfy this ordering is 



X)s=o (s) = 2^- Let us assume that the Xj, too, are pairwise distinct. Then we 
can order them as Ai < . . . < Ap, because of (|95|) . The total number of states 
satisfying the above restrictions is J2p=(} (p)^^ — S'^, which is the dimension 
of our Hilbert space. This strongly suggests that our algebraic Bethe ansatz 
solution provides us with a complete set of states. An actual proof would 
require to calculate all scalar products of the form (Ai, . . . , Xp\X[, . . . , A^), 
which is feasible, yet beyond the scope of this article. 

We want to point out, that the states |Ai, . . . , Ap) are eigenstates of the 
z-component of the total spin operator. The proof goes as follows. The R- 
matrix of the Hubbard model in the infinite coupling limit is invariant under 
Co — Co). Now 



2 



^^ = iEKT-^.i) = iE(^.2-^.3) ' (123) 

j=l j=l 

and eq. (|iB| ) implies that [S^,T{u)] = ^[T{u),e1 — Cg], or, for the elements 
Ca{u) of the monodromy matrix, 

[S^CM] = l{aXaiu) . (124) 
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It follows that 

S^\X^,...,X,) = i(S^Fr-'^-C„,(Ai)...a^(Ap)|0) 

= i(p-2s)|Ai,...,Ap) . (125) 

Let us compare our results with the recent coordinate Bethe ansatz so- 
lution of the model by Izergin et al. |22|. The one-particle algebraic Bethe 
ansatz states are easily shown to be of the form 

C,(A)|0)ocX:(-tan(A)y4jO) , (126) 
i=i 

where Cja = Cj^, if a = 1, and Cja = Cji, if a = 2. From this equation we 
identify the quasimomenta k{X) of one-particle states as 

ei'^W = -tan(A) . (127) 

Let us perform a similar reparametrization for the auxiliary spin problem, 

ei'?('^) = tan(^) . (128) 

Then, using ( p.l5| ), we can rewrite the Bethe ansatz equations (|102|) and ( |lllj ) 

as 

gifc.i = giELi^fe , j = l,...,p , (129) 
gigfeP = ^ k=l,...,s , (130) 

where kj = k{Xj), qu = <i{l^k)- These equations agree with the Bethe ansatz 
equations obtained in p2| . 

Eq. ( pTl) implies that the energy of a state |Ai, . . . Ap) is 

E = 5„ln(A(w))|„^o = X:(tan(A,)+cot(A,)) 

i=i 

= -2X:cos(fc.) . (131) 

Eqs. ( |122| ) and ( |125|) show that the states |Ai, . . . Xp) are also eigenstates of 
the particle number operator and the ^-component of the total spin S^. We 
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may thus add a chemical potential and a magnetic field to our Hamiltonian 

L 

Hf,,B = ^{-A(c]^Cj+i,^+c]+i^^Cj^)A - /iKr+nji) + B{nj^-nji)} . (132) 
i=i 

Then if^,B|Ai, . . . , Ap) = -E;,,b|Ai, . . . , Ap), where 

p 

= -2^cos(A;,) -/ip + 5(p-2s) . (133) 
j=i 

This means that the algebraic Bethe ansatz states |Ai, . . . , Ap) can be used 
to calculate form factors and correlation functions in the grand canonical 
ensemble (cf @). 



Yangian symmetry of the Hubbard model in the infinite 
coupling limit 

The Hubbard Hamiltonian on the infinite interval is invariant under the direct 
sum of two su(2) Yangian quantum groups p3|, We would like to point 
out here, that one of these Yangian symmetries survives the infinite coupling 
limit. First of all, let us recall |2^ that the su(2) Yangian is a Hopf algebra 
which is spanned by two triples of generators 1°^, J", a = x,y,z, satisfying 
the relations 

[r,l'] = CabJ' , (134) 
[r,j'] = CabcJ' , (135) 

[[r,j\[r,.f]] + [[J^J1,[/^J^]] = 

-^{flabefghCcde + a'cdefghCabe){I^ , , J^} ■ (136) 

Here Cabc = i^"'^'^ is the antisymmetric tensor of structure constants of su(2), 
and aabcdef = CadgCbehCcfiCghi- The bracket { } in ([1361) denotes the sym- 
metrized product 

{xi,X2,X3} = ^ ^ XiXjXk . (137) 



30 



Being a Hopf algebra Y(su(2)) carries an outer structure (co-multiplication, 
antipode, co-unit), which guarantees that Y(su(2)) has a rich representation 
theory 



(138) 
(139) 



As a corollary of a theorem proven in '^Tj it follows that 



form a representation of Y(su(2)). The Sj in this equation are spin density 
operators in fermionic representation (|H^). obviously commutes with the 
Hubbard Hamiltonian in the infinite coupling limit (pOf). We will give a simple 
proof I^H] that J'*, too, commutes with the Hamiltonian (pOD, if we replace the 



summation in (0) by a summation over all integers. Let Sjj^ 



2 a, 



and define 



It was shown in j^S], that the Hubbard Hamiltonian Hh, 



Hh 
T 



T + UD 

-E(4 



(140) 

(141) 
(142) 



(143) 



commutes with 



ja 



(144) 



and that J" and form a representation of the Y(su(2)) Yangian. Since 
commutes with Hh for all real U, we obtain the identity [T, J"] = [K^-jD]. 
Now DA = AD = 0, and [A, J""] = 0, since [A, S^] = 0. It follows that 

[H, r] = [ATA, r] = A[T, J'^]A = A[K% D]A = . (145) 

Let us remark that a more systematic way to explore the Yangian sym- 
metry of the model would be to apply the approach developed in |2^, ^], 
where the Yang-Baxter algebra was investigated in the thermodynamic limit. 
This approach would also provide us with the action of the Yangian on eigen- 
states. 
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Conclusions 

We would like to stress that our approach is widely applicable. Given a 
solution R{u, v) of the Yang-Baxter equation (^2]), it consists of the following 
steps. 



(i) Use eq. (|26|) to choose a grading which is compatible with the i?-matrix. 
In general, this choice is not unique. A given i?-matrix may be com- 
patible with different gradings. 

(ii) For a given grading choose a fermionic representation. Start with a 
tensor product representation (^) of projection operators Yj^ of suf- 
ficiently high dimension and adjust it to the grading by deleting rows 
and columns of 1^-^. Again the choice is not unique. 



(iii) Write down the Hamiltonian (^OD and replace the graded spin operators 
^ja fermionic projection operators Yj^. This Hamiltonian is inte- 
grable by construction. There may be different methods to diagonalize 
it. 

(iv) The Hamiltonian can be diagonalized by (nested) coordinate Bethe 
ansatz. 

(v) It may be possible to diagonalize the Hamiltonian algebraically. In 
this case the L-matrix is obtained from eq. ([27|), and the graded Yang- 



Baxter algebra is given by (|30|). Note however, that there is no general 
recipe for an algebraic Bethe ansatz. Depending on the structure of 
the i?-matrix and of the monodromy matrix an algebraic Bethe ansatz 
may be a difficult task (cf e.g. [^, 0), if it is possible at all. 



Our approach circumvents certain difficulties connected with the Jordan- 
Wigner transformation, which, in certain cases, may be alternatively used to 
connect a spin representation of an integrable system with a fermionic repre- 
sentation. In our approach there is no twisting of the boundary conditions, 
and symmetries are directly obtainable via eq. (H). Our approach works for 



higher su(iV) spins where no Jordan- Wigner transformation is known. 

The main example to illustrate our ideas was the Hubbard model in the 
infinite coupling limit. An interesting lesson to learn from this example is 
that the Hubbard model in the infinite coupling limit is equivalent to a su(3)- 
spin generalization of the XX spin chain. 
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